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Abstract 

A  technique  has  been  developed  for  the  calculation  of  excited  state, 
one-electron  wave  functions  based  on  the  Thomas-Fermi  statistical  theory  of 
the  atom.  The  technique  is  applicable  to  heavy  atoms  for  vhich  Hartree  type 
solutions  are  complex  and  difficult  to  obtain.   In  this  paper  the  previously 
obtained  Thomas-Fermi  core  potential  for  the  cesium  atom  and  a  Heisenberg  type 
polarization  correction  is  used  as  a  central  field  in  the  Schrodinger  equation. 
Correction  for  penetration  of  the  excited  electron's  orbital  is  made,  and  the 
Biermann-Lubeck  approach  for  solving  the  wave  equation  is  utilized.  This  allows 
for  the  inclusion  of  a  qualitative  correction  for  exchange.  The  cesiiun  atom's 
7s  excited  state  which  has  not  been  obtained  by  any  Hartree  method  is  computed. 
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1.   Introduction. 

The  only  systematic  approach  to  the  calculation  of  atomic  structure  and.  prop- 
erties is  via  the  Hartree  ^ J  or  Hartree-Fock  L J  approximation.   Hovever^  there 
are  limitations  to  these  calculations .  Their  numerical  complexity  makes  it  diffi- 
cult to  obtain  an  overall  physical  picture  of  the  atom.   Only  isolated  solutions 
are  possible  in  these  approximations  and  for  heavy  atoms  the  numerical  complexity 
of  the  self  consistent  iterations  makes  the  calculations  for  atomic  structure  all 
but  impossible. 

There  are  additional  complications  which  arise  when  attempting  to  calculate 
the  structure  of  excited  states  since  one  requires  that  the  wave  functions  used  be 
orthogonal  to  all  lower  states^  and  it  is  sometimes  Impossible  to  obtain  solutions 
with  this  requirement. 

For  those  states  for  which  Hartree  or  Hartree-Fock  one-electron  wave  functions 
cannot  be  obtained^  one  can  use  a  somewhat  cruder  approximation  which  has  the  ad- 
vantage  of  simplicity,  namely  the  Thomas-Fermi  modell-'— J.   The  basic  assumptions 
In  this  model,  however,  require  that  we  investigate  only  the  ground  state  conflgior- 
ations.   Some  attempts  have  been  made  by  Latter  L'-J  to  investigate  the  excited  states 
of  an  electron  in  the  Thomas-Fermi  central  potential.  This  is  of  course  an  approx- 
imation to  the  description  of  excited  states  of  atoms  since,  for  example,  some 
self -force  on  the  electron  is  included  and  there  is  no  correction  for  interaction 
effects.  A  procedure  has  been  previously  described  by  this  author L ^ U  vhlch  Is 
a  combination  of  the  Thomas-Fermi  technique  and  a  self-consistent  procedure.   It 
is  easiest  to  apply  in  cases  where  there  are  several  closed  shells  and  the  excita- 
tions being  discussed  are  those  of  electrons  outside  these  closed  shells.   The  ex- 
cited electron  is  considered  to  move  in  the  Thomas-Fermi  field  of  the  remaining 
electrons,  and  a  self-consistent  procedure  is  used  to  account  for  the  distortion 
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of  the  Thomas -Fermi  field  due  to  the  penetration  of  the  charge  cloud  by  the  ex- 
cited electron.   In  this  report,  the  solution  to  the  7s  excited  state  of  the 
cesium  atom  is  presented. 
2.   The  General  Technique. 

For  an  alkali  atom,  we  first  neglect  the  interaction  of  the  excited  electron 
vith  the  remaining  ion  core.  The  ion  core  can  then  be  treated  by  the  Thomas-Fermi 
technique  since  the  remaining  core  electrons  are  all  in  the  ground  state^  and  for 
a  heavy  ion,  well  approximated  by  a  continuous  distribution.   Such  an  ionic  system 
obeys  the  following  differential  system  ^  -■ 

with   ^{0)   =  1,  (2.2) 

<^ixj   =  0,  (2.3) 

X  ?^'(x^)  =  -(Z-N)/Z  (2.J+) 

where  <^{r)   =  r rv(r)-V  ~j/Ze  (2.5) 

X  =  r/a^(9rt^/l28z)^/5  =  rZ^/^/0.8855i^la^  (2.6) 

with  ■V(r),  the  radially  symmetric  potential;  a  ,   the  first  Bohr  radius;  Z,  the 

nuclear  charge;  N,  the  number  of  electrons  in  the  ion  core;  and  r  ,  the  effective 

radius  of  the  ion,  being  defined  such  that  for  r>r  ,  ■V(r)  =  (Z-N)e/r. 

Solving  this  set  of  equations  yields  a  potential  distribution  for  the  core 

which  will  be  called  V  .       .  The  solution  of  the  Thomas-Fermi  ion  has  been 

ion  core 

given  in  reference  [6] . 

The  potential  will  not  be  exactly  the  ion  that  the  excited  electron  moves 
in.  'Fr-ren   for  a  completely  non- penetrating  outer  electron,  which  might  be  expected 
to  move  in  a  completely  hydrogen-like  Coulomb  field,  there  is  produced  a  resultant 
asymmetry  or  polarization  which  leads  to  slightly  greater  binding  energy  for  the 
nonpenetrating  electron  than  for  the  hydrogenic  electron.   In  the  present  analysis 


r9i 

this  correction  is  made  using  a  term  introduced  by  Born  and  Heisenberg  •- -" ,  the 
argument  being  presented  in  reference  [7] .  The  latter  shows  there  exists  a  polar- 
ization energy  of 

?   I4- 
U  ,  =  -a  e  /2r  ,  (2.?) 

pol     P  \   '  / 

where  a  is  the  polarizability  of  the  ion  core.   Its  value  for  the  cesium  ion  is 

-2k 
approximately  2.42  in  10    units. 

3.   Solution  of  the  Schr6'dinger  Equation. 

The  charge  distribution  and  consequently  the  potential  of  the  ion  can  now  be 
used  as  a  central  field  in  which  the  valence  electron  moves.  The  calculated  ex- 
cited and/or  ground  state  one-electron  wave  functions  will  give  the  charge  dis- 
tribution of  the  valence  electron  for  the  states  considered.   This  statement  is 
only  approximately  true  for  several  reasons,  the  most  important  of  which  is  that 
we  have  calculated  a  charge  distribution  for  the  core  which  is  good  only  in  the 
Thomas-Fermi  approximation.   Other  approximations  which  will  have  to  be  made 
concern  the  partial  shielding  of  the  core  due  to  the  penetration  of  the  valence 
electron's  charge  distribution,  the  exchange  effects,  and  the  dynamical  aspects 
of  the  problem  which  include  the  deformation  of  the  core  due  to  the  presence  of 
the  penetrating  electron.   It  will  be  seen  that  these  last  three  effects  can  be 
included  approximately  in  our  calculation. 

If  we  consider  a  separable,  product  type  wave  function,  and  make  the  further 
assimption  that  angular  dependence  can  be  separated  out,  the  one-electron  wave 
function  of  the  1   electron  is  given  by 

*1  =  K(^)Yt,Je'5^)  =  P'(^)A]\ni(e,?()  (3.1) 

and  the  wave  equation  for  the  valence  electron  moving  in  the  field  of  the  Thomas- 
Fermi  potential  of  the  ion  is 


x;'(r) 
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(2me./ft2)  -  (2meV^(r)/fi2)  +  (t)(t+l)/r^ 


X^(r)  =  0       (5.2) 


where  Y.   (6,^)  are  the  spherical  harmonic  fimctions,  V   (r)  the  core  potential^ 
I   the  angular  momentimi  quantum  number,  and  e .  the  corresponding  eigenvalue .   In 
accordance  with' normalization  we  further  require  that 


rr^(r)R*(r)dr  =  j     X(r)X*  (r)dr  =  1.  (5.5) 

o  o 

It  is  seen  from  the  above  that  after  angular  dependence  is  separated  out, 
the  remaining  second-order  equation  can  be  put  into  a  fonn  that  does  not  contain 
the  first  derivative  by  introducing  the  X(r)  function.   It  is  then  amenable  to 
numerical  integration  by  the  Milne '-   -I  method  or  by  standard  numerical  techniques 
due  to  HartreeL  J,  Gauss-Jackson-Jumerov'—  '    -^ ,    and  Blanch*-  J. 

As  indicated  in  Section  2,  the  Heisenberg  form  for  the  polarization  potential 

2/  h 
energy  is  -a  e  /2r  ,  where  the  argument  for  its  determination  holds  only  in  the 


region  outside  the  core.   If  this  term  were  added  only  in  the  outer  region,  there 


would  be  a  discontinuity  at  the  core  boundary.  Following  Biermann  and  Hartlng 

2 


for  the  inner 


we  add  a  polarization  term  of  the  form  (a  e  /2r  ) 

region  as  well.   Here  f  is  a  fraction  which  was  taken  as  0.^  for  this  calculation. 

Although  the  exponent  can  be  raised  to  any  power  greater  than  foiir,  eight  was 

ri5i 

found  to  give  the  best  fit  with  optical  values  "^  <J  .   The  addition  of  this  term 

has  the  effect  of  eliminating  the  polarization  term  in  the  core  until  r  =  fr  . 

/  h 
It  then  goes  smoothly  into  the  1/r  behavior  at  the  core  boundary.  Although 

strictly  speaking  this  is  an  empirical  correction  it  is  a  small  one,  and  too 

much  concern  need  not  be  given  to  the  exact  form  of  the  correction. 

A  useful  apporach  to  the  solution  of  the  wave  equation  is  that  developed 

by  Biermann  and  Lubeck ' '.   In  addition  to  the  polarization  contribution,  and 
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being  adaptable  to  the  Thomas-Fermi  potential,  the  method  also  allows  for  the 
inclusion  of  a  semi-empirical  exchange  correction.   Introduction  of  the  latter 
correction  has  the  effect  of  introducing  a  non- experimentally  determined  para- 
meter into  a  system  which  is  already  determined  without  redundancy.   In  order 
to  compensate  for  this  while  calculating  the  best  possible  wave  function  for 
the  assumed  potential  (originally  taken  as  the  Hartree  potential),  Biermann  and 
Liibeck  solve  the  wave  equation  utilizing  the  experimentally  determined  term  value, 
Accurate  listings  of  these  values  can  presently  be  found  in  Circular  h-67   of  the 
National  Bureau  of  Standards  L  J .   Starting  with  this  experimental  term  value  is 
not  necessarily  an  unfair  approach  since  both  Ifertree  and  Thomas-Fermi  potentials 
do  lead  to  term  values  which  are  very  close  to  the  experimentally  determined  ones, 

The  exchange  correction  made  by  Biermann  and  Liibeck  can  be  applied  here  by 
multiplying  the  Thomas-Fermi  central  potential  distribution  by  the  function 

B(r)  =  (1  +  3re"^''/^''o^  )   ;  (5-4) 

for  the  computation  it  is  convenient  to  express  r  and  P  coefficients  of  the  ex- 

o       o 
ponential  in  A  and  1/A  units,  respectively.  With  this  approach  the  radial  wave 

equation  to  be  solved  is  now 


X'  '  - 


2me 


a  2 


ex£  _  2mgV(r)  ^  (t)(t+l)  _   pe^   ^^_^-(r/fr^)  ^ 


^2 


where  V(r)  =  rZesz((r)/r  -  e/r~\ 


2r 


,  „  -(r/fr  y 
1+Pre  ^    '      o' 


X  =  0 


=  0  ;   0  S  r  g  r 


(3.5) 


(3.6) 


and    V(r)  =  (e/r) 


l+Pre-(^/^^o^ 


^ 


r  >  r 


(3.7) 


The  solution  of  the  equation  for  large  values  of  r  is  given  as  X(r)  =  ke 


±\|2me.r/^. 


Varying  p  for  a  particular  eigenvalue  and  angular  momentum  state  leads  to  solutions 
such  that  all  diverge  as   e    except  the  one  associated  with  the  desired  value  of 
p;  in  this  case  we  have  the  physically  required  exponential  decay. 

Since  this  procedure  takes  advantage  of  the  experimentally  determined  term 
value^  and  makes  at  least  partial  correction  for  almost  all  other  effects,  it  can 
be  expected  to  yield  fairly  accurate  wave  functions,  even  for  excited  states.   In- 
deed this  is  the  case  and  Biermann  and  Lttbeck  find  that  oscillator  strengths  for 
the  states  in  the  principal  series  agree  very  well  with  those  obtained  from  experi- 
mental measurements.  The  Biermann  and  L^beck  procedijire  neglects,  however,  the  ef- 
fect of  the  penetration  of  the  excited  electron  into  the  core,  and  the  modification 
of  the  core  potential  due  to  this  penetration.   While  this  is  a  small  correction 
in  the  case  of  the  more  highly  excited  levels,  it  can  be  significant  for  the  first 
few  excited  states.   It  is  possible  to  include  the  effect  of  this  penetration  by- 
utilizing  a  self  consistent  iterative  process  described  in  reference  [7]^  and  sum- 
marized in  the  next  section. 
h.      The  Self  Consistent  Technique. 

The  effect  of  the  penetration  of  the  valence  electron  into  the  core  is  parti- 
ally to  shield  the  nuclear  charge  and  consequently  to  increase  the  size  of  the  ion 
core.   The  form  of  the  core  charge  distribution  also  changes,  and  with  it,  the  ef- 
fective potential  in  which  the  valence  electron  moves  changes.   Since  the  effective 
potential  depends  on  the  charge  distribution,  and  this  in  turn  depends  on  the  ef- 
fective potential,  we  have  in  principle  the  ingredients  of  a  self  consistent  pro- 
cedure. 

We  first  solve  the  problem  neglecting  penetration  of  the  core.   If  the  portion 
of  the  computed  esLclted  electron's  charge  distribution  penetrating  the  core  Is  a. 
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the  core  size  increases  somewhat,  and  the  potential  at  the  edge  of  the  core  be- 
comes  (Z-N-a)e/r  ,   where  r'  Is  the  new  core  radius.  This  new  core  potential 
combined  with  the  polarization  contribution  is  used  to  solve  the  Schrodinger  eq- 
uation for  the  new  charge  distribution  until  the  solution  is  self  consistent. 

There  are  complications  in  carrying  out  this  procedure.  The  penetration 
not  only  leads  to  a  modified  set  of  boundary  conditions,  but  strictly  speaking, 
to  a  modified  Thomas-Fermi  equation  as  well. 

The  Thomas-'Fermi  potential  distribution  for  the  pure  ion  represents  a  zeroth 

order  approximation  to  the  ion  core.  We  can  use  this  potential  with  a  polarization 

correction  in  the  Schrodinger  equation  to  calculate  an  energy  eigenvalue  and  a 

wave  function  for  the  excited  electron  e  which  is  in  some  state  characterized 

by  specific  principal  and  azimuthal  quantum  niunbers.  From  the  one  electron  wave 

function,  ijj  ^,  so  determined,  we  can  calculate  a  number  density  for  the  excited 

e 
electron 


n  * 


t 


e 


2 


(^.1) 


We  can  also  define  a  dimensionless  potential,  that  is  calculated  from  the  potential 

2      , 
V  ^  as  determined  by  Poisson's  equation,  ^Y  ^  =  4jt  n  ^e: 

e  e        e 

9^  ^  =  rV  ^(r)/Ze.  {k .2) 

e      e 

It  is  shown  in  reference  [tI ,  that  if  this  potential  is  taken  into  account,  the  core 
potential,  </>rr^,   must  now  satisfy  a  modified  Thomas-Fermi  equation, 

e 


As  in  the  case  for  the  pure  ion,  the  boxindary  conditions  at  the  nucleus  and  limit- 


9  - 


ing  radius  are 

5^^(x)  =1  ;   ^^{xj   =  0.  (4.M 

The  detail,  but  not  the  physical  content,  of  the  third  boundary  relation  alters. 
It  becomes 

-(Z-N-a)/z  =  ^J^{^^)      ■  (^-5) 

5.   Starting  Values. 

The  inner  region  starting  values  can  be  obtained  by  expanding  X(r)  as  a 
power  series  in  r,  substituting  into  the  differential  equation  and  setting  the 
collected  coefficients  of  the  r  powers  to  zero.   If  exchange  effect  near  the  origin 
is  neglected,  we  can  write  the  Schrodinger  equation  as 


a   +   s^/r  -  s^<^{r)/v 


X(r)  =  0  (5.1) 


X"(r)  - 

where 

CT  =  2m/-h^(€-eV  )   ,  s^  =  (e)(l+l)    ,  and  s  =  2me  Z/ft^.  (5-2) 

Expressing  the  wave  function  in  power  series  form 

CO 

X(r)  =  )  a^rn    ,  (5.5) 

n=0 

and  the  collected  coefficient  of  the  general  term  r   is 

(m+l)(m+2)a^_2  -  oa^   +  S2i/(r)a^_^  -  s^a^_2  =  0  .  (5.^ 

The  5z((r)  function  varies  between  zero  and  unity  through  the  entire  inner  region 
and  is  treated  as  a  parameter  in  the  expansion.  The  series  must  be  evaluated  from 
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m  =  -2  since  there  is  an  inverse  square  term.   Combining  the  above  results  for 
the  various  angular  momenta  yields  the  following  expansion  for  the  s  state: 


X        (r) 

3- state  ^  ' 


r  -  (S2^(r)/2)r^  +  (l/6) 


(a  +  bI  f   (r)/2^ 


■(l/l2)(sJ(r)/2) 


a  +  (l/3)(a  +  si   9^2(r)/2) 


+(1/20)  J  {a/G){a   +  sI<^^{t)/2) 


slf{T)/{2k) 


a  -   (l/5)(a  +  sl^fyr)/^    I  r^ 


(5.5) 


Even  neglecting  the  exchange  correction  term,  the  coefficients  are  seen  to 
become  complex  as  higher  powers  are  needed.   Unfortunately  It  will  be  seen  that 
heavy  atom  wave  functions  are  rapidly  changing  in  the  strong  Coulomb  field  near 
the  nucleus^  and  even  if  a  fairly  small  interval  is  taken  for  the  X(r)  tabulation, 
using  the  Milne  method  of  integration,  the  starting  Amlue  calculation  becomes 
tedious.  Although  the  present  calculation  was  carried  out  in  this  manner,  the 
authors  recommend  the  Gauss-Jackeon-Numerov'-  '    -^   method  for  future  work.   In 
that  case,  only  two  or  three  terms  in  the  expansion  in  (5-5)  are  needed.   The 
approach  requires  that  the  Thomas-Fermi  potential  be  determined  at  Irregular 
intervals  near  the  origin,  but  that  is  less  cumbersome  than  determining  the  start- 
ing values  using  a  fifth  or  sixth  degree  polynomial. 

It  can  also  now  be  pointed  out  that  the  parameter  in  the  Biermann  and  Lubeck 
scheme  must  be  known  with  ever  greater  accirracy  as  heavier  atoms  are  studied.   In 
their  original  work,  Biermann  and  Lubeck  required  only  three  significant  figures 
for  convergence.   D.  Villars L  -■  noted  that  five  figures  were  required  for  his  work 
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with  potassium.   Unfortimately,  for  the  55  electron  cesium  atom,  eight  significant 
figures  were  foiind  to  he  necessairy.  This  begins  to  exhaust  the  capabilities  of 
even  advanced  computers,  and  the  author  suggests  that  for  heavier  atoms,  the 

ri9i 

Ridley  I-  -'  scheme  be  used  since  the  exchange  correction  is  smaller. 

6.  Results. 

As  discussed,  the  7s  state  of  the  relatively  heavy  cesium  atom  was  calculated 
using  the  above  procedure.  This  state  is  only  one  excited  level  above  the  ground 

state  and  therefore  shows  a  fair  degree  of  penetration.   Using  the  Biermann  and 

-12 
Lubeck  computing  approach,  with  the  experimental  term  value  of  2.556^39  x  10 

ergs,  a  final  degree  of  penetration,  a,  was  determined  as  'J.O'ji,.     Figure  1  pre- 
sents the  determined  wave  function.  Figure  2  shows  the  inner  region  expanded, 
and  Figure  5  illustrates  the  probability  distribution.  Table  1  presents  a 
summary,  and  Table  2  presents  the  full  detail  of  the  calculation. 
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EPSILON 


2.556439E-12 


RMAX 


2.450235E-07 


INNER  SOLUTION 


CHI 


NORMALIZATION 
8.246397E  00 

CHISQ 


2.357331E- 

-11 

1.515622E 

02 

2.297110E 

04 

1.023573E- 

-09 

-9.890845E 

02 

9.782882E 

05 

2.023573E- 

-09 

1.031266E 

03 

1.063510E 

06 

3.023573E- 

-09 

4.474262E 

02 

2.001902E 

05 

4.023573E- 

-09 

-8.278346E 

02 

6.853101E 

05 

5.023573E- 

-09 

-1.337223E 

03 

1.788165E 

06 

6.023573E- 

-09 

-1.038766E 

03 

1.079035E 

06 

7.023573E- 

-09 

-3.096236E 

02 

9.586b80E 

04 

8.023573E- 

-09 

5.037750E 

02 

2.537892E 

05 

9.023573E- 

-09 

1.186158E 

03 

1.406972E 

06 

1.002357E- 

-08 

1.637021E 

03 

2.679838E 

06 

1.102357E- 

-08 

1.832169E 

03 

3.356843E 

06 

1.202337E- 

-08 

1.791166E 

03 

3.208274E 

06 

1.302357E- 

-08 

1.555250E 

03 

2.418804E 

06 

1.402357E- 

-08 

1.173848E 

03 

1.377919E 

06 

1.502357E- 

-08 

6.968905E 

02 

4.856564E 

05 

1.602357E- 

-08 

1.707055E 

02 

2.914035E 

04 

1.702357E- 

-08 

-3.640405E 

02 

1.325255E 

05 

1.802357E- 

-08 

-8.732911E 

02 

7.626373E 

05 

1.902357E- 

-08 

-1.329982E 

03 

1.768853E 

06 

2.002357E- 

-08 

-1.714144E 

03 

2.938290E 

06 

2.102357E- 

-08 

-2.012828E 

03 

4.051478E 

06 

2.202357E- 

-08 

-2.219864E 

03 

4.927795E 

06 

2.302357E- 

-08 

-2.335359E 

03 

5.453901E 

06 

2.402357E- 

-08 

-2.364870E 

03 

5.592611E 

06 

2.502357E- 

-08 

-2.318172E 

03 

5.373922E 

06 

2.602357E- 

-08 

-2.207670E 

03 

4.873806E 

06 

2.702357E- 

-08 

-2.046629E 

03 

4.188688E 

06 

2.802357E- 

-08 

-1  .847520E 

03 

3.413331E 

06 

2.902357E- 

-08 

-1 .620835E 

03 

2.627107E 

06 

3.002357E- 

-08 

-1.374587E 

03 

1.889489E 

06 

3.102357E- 

-08 

-1.114507E 

03 

1.242126E 

06 

3.202357E- 

-08 

-8.446624E 

02 

7.134545E 

05 

ALPHA(RO) 
7.070833E-02 
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OUTER  SOLUTION 

R 

CHI 

CHISO 

3.242352E-08 

-7.346993E 

02 

5.397830E 

05 

4.2^2352E-08 

1.975877E 

03 

3.904092E 

06 

5.242352E-08 

3.828054E 

03 

1.465400E 

07 

6.242352E-08 

4.544958E 

03 

2.065664E 

07 

7.242352E-08 

4.409886E 

03 

1.944709E 

07 

8.2A2352E-08 

3.808746E 

03 

1.450655E 

07 

9.242352E-08 

3.044773E 

03 

9.270641E 

06 

1.024233E-07 

2.302172E 

03 

5.299995E 

06 

1.124235E-07 

1.668634E 

03 

2.784339E 

06 

1.22A235E-07 

1.169875E 

03 

1.368607E 

06 

1.324235E-07 

7.984085E 

02 

6.374561E 

05 

1.A24235E-07 

5.329000E 

02 

2.839824E 

05 

1.524235E-07 

3.490897E 

02 

1.218636E 

05 

1.624235E-07 

2.250614E 

02 

5.065263E 

04 

1.724235E-07 

1.431219E 

02 

2.048388E 

04 

1.82A235E-07 

8.994511E 

01 

8.090123E 

03 

1.924235E-07 

5.596465E 

01 

3.132042E 

03 

2.024235E-07 

3.455357E 

01 

1.193949E 

03 

2.124235E-07 

2.125141E 

01 

4.516222E 

02 

2.224235E-07 

1.313141E 

01 

1.724340E 

02 

2.324235E-07 

8.323815E 

00 

6.928589E 

01 

2.424235E-07 

5.680861E 

00 

3.227218E 

01 
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EPSILON  RMAX 

2.556439E-12        2.450235E-07 
INNER  SOLUTION 


CHI 


2.357331E-11 
1.235733E-10 
2.235733E-10 
3.235733E-10 
4.235733E-10 
5.235733E-10 
6.235733E-10 
7.235733E-10 
8.235733E-10 
9.235733E-10 
1.023573E-09 
1.123573E-09 
1.223573E-09 
1.323373E-09 
1.423573E-09 
1.523573E-09 
1.623573E-09 
1.723573E-09 
1.823573E-09 
1.923573E-09 
2.023573E-09 
2.123573E-09 
2.223573E-09 
2.323573E-09 
2.423573E-09 
2.523573E-09 
2.623573E-09 
2.723573E-09 
2.823573E-09 
2.923573E-09 
3.023373E-09 
3.123573E-09 
3.223573E-09 
3.323573E-09 
3.423573E-09 
3.523573E-09 
3.623573E-09 
3.723573E-09 
3.823573E-09 
3.923573E-09 
4.023573E-09 
4.123573E-09 
4.223573E-09 
4.323373E-09 
4.423573E-09 
4.523573E-09 
4.623572E-09 
4,723372E-09 
4.823572E-09 
4.923572E-09 
5.023572E-09 
5.123572E-09 
5.223572E-09 
5.323572E-09 


NORMALIZATION 
8.246397E  00 

CHISQ 


ALPHA(RO) 
7.070833E-02 


1 

.515622E 

02 

2.297110E 

04 

1 

.990095E 

02 

3.960477E 

04 

■1 

.011668E 

02 

1.023472E 

04 

-3 

.479303E 

02 

1.210555E 

05 

■5 

.931107E 

02 

3.517802E 

05 

•8 

.240385E 

02 

6.790394E 

05 

•1 

.019018E 

03 

1.0383^8E 

06 

•1 

.146263E 

03 

1.313918E 

06 

■1 

,184067E 

03 

1.402015E 

06 

■1 

.128031E 

03 

1.272453E 

06 

•9 

.890845E 

02 

9.782882E 

05 

•7 

.870408E 

02 

6.194332E 

05 

■5 

.444819E 

02 

2.964605E 

05 

•2 

.827089E 

02 

7.992432E 

04 

•1 

.973898E 

01 

3.896271E 

02 

2 

.303196E 

02 

5.304710E 

04 

4 

.571483E 

02 

2.089845E 

05 

6 

.537706E 

02 

4.274160E 

05 

8 

.160009E 

02 

6.658575E 

05 

9. 

.418975E 

02 

8.871709E 

05 

I. 

.031266E 

03 

1.063510E 

06 

1. 

.085231E 

03 

1.177726E 

06 

1. 

.105869E 

03 

1.222947E 

06 

1. 

.095917E 

03 

1.201034E 

06 

1. 

.058528E 

03 

1.120482E 

06 

9. 

.970825E 

02 

9.941735E 

05 

9. 

.150408E 

02 

8.372997E 

05 

8. 

.158314E 

02 

6.655809E 

05 

7. 

027680E 

02 

4.938829E 

05 

5. 

.789908E 

02 

3.352303E 

05 

4. 

.474262E 

02 

2.001902E 

05 

3. 

107624E 

02 

9.657327E 

04 

1. 

714361E 

02 

2.939033E 

04 

3. 

162829E 

01 

1.000349E 

03 

1. 

067328E 

02 

1.139189E 

04 

2. 

419641E 

02 

5.854661E 

04 

3. 

726164E 

02 

1.388430E 

05 

4i 

974607E 

02 

2.474671E 

05 

6. 

154730E 

02 

3.788071E 

05 

7. 

258184E 

02 

5.268124E 

05 

8. 

278346E 

02 

6.853101E 

05 

9. 

210159E 

02 

8.482703E 

05 

004998E 

03 

1.010020E 

06 

079541E 

03 

1.165408E 

06 

144518E 

03 

1.309921E 

06 

199900E 

03 

1.439759E 

06 

245742E 

03 

1.551874E 

06 

282176E 

03 

1.643976E 

06 

309395E 

03 

1.714516E 

06 

327646E 

03 

1.762645E 

06 

337223E 

03 

1.788165E 

06 

338454E 

03 

1.791459E 

06 

331699E 

03 

1.773423E 

06 

317342E 

03 

1.735390E 

06 
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5.423 
5.523 
5.623 
5.723 
5.823 
5.923 
6.023 
6.123 
6.223 
6.323 
6.423 
6.523 
6.623 
6.723 
6.823 
6.923 
7.023 
7.123 
7.223 
7.323 
7.423 
7.523 
7.623 
7.723 
7.823 
7.923 
8.023 
8.123 
8.223 
8.323 
8.423 
8.523 
8.623 
8.723 
8.823 
8.923 
9.023 
9.123 
9.223 
9.323 
9.423 
9.523 
9.623 
9.723 
9.823 
9.923 
1.002 
1.012 
1.022 
.032 
.042 
.052 
.062 
.072 
.082 
.092 
.102 
.112 


1.122 
1.132 
1.142 
1.152 


572E-09 
572E-09 
572E-09 
572E-09 
572E-09 
572E-09 
572E-09 
572E-09 
572E-09 
572E-09 
572E-09 
572E-09 
572E-09 
571E-09 
571E-09 
571E-09 
571E-09 
571E-09 
571E-09 
571E-09 
571E-09 
571E-09 
571E-09 
571E-09 
571E-09 
571E-09 
571E-09 
570E-09 
570E-09 
570E-09 
570E-09 
570E-09 
570E-09 
570E-09 
570E-09 
570E-09 
569E-09 
569E-09 
569E-09 
569E-09 
569E-09 
569E-09 
569E-09 
569E-09 
569E-09 
569E-09 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 
357E-08 


■1.295783E    03 

■1.267437E    03 

■1.232726E    03 

•1.192076E    03 

■1.145916E    03 

■1.094672E    03 

■1.038766E    03 

■9.786137E    02 

-9.146221E    02 

•8.471888E    02 

-7.767000E    02 

■7.035300E    02 

-6.280401E    02 

•5.505779E    02 

■4.714773E    02 

•3.910576E    02 

•3.096236E    02 

-2.274656E  02 

-1.448590E  02 

■6.206492E    01 

2.067004E  01 

1.031134E  02 

1.850466E  02 

2.662646E  02 

3.465759E  02 

4.258014E  02 

5.037750E  02 

5.803423E  02 

6.553609E  02 

7.286993E  02 

8.002372E  02 

8.698643E  02 

9.374806E  02 

1.002995E  03 

1.066327E  03 

1.127403E  03 

1.186158E  03 

1.242537E  03 

1.296488E  03 

1«347972E  03 

1.396951E  03 

1.443397E  03 

1.487288E  03 

1.528605E  03 

1.567336E  03 

1.603476E  03 

1.637021E  03 

1.667974E  03 

1.696342E  03 

1.722135E  03 

1.745366E  03 

1.766055E  03 

1.784220E  03 

1.799886E  03 

1.813080E  03 

1.823831E  03 

1.832169E  03 

1.838129E  03 

1.841746E  03 

1.843058E  03 

1.842103E  03 

1.838924E  03 


.679055E 
.606397E 
.519613E 
.421045E 
.313123E 
.198307E 
.079035E 
.576848E 
.365337E 
.177289E 
.032629E 
.949545E 
.944343E 
.031361E 
.222909E 
.529261E 
.586680E 
.174059E 
.098413E 
.852054E 
.272504E 
.063237E 
.424224E 
.089686E 
.201148E 
.813069E 
.537892E 
.367972E 
.294979E 
.310027E 
.403796E 
.566640E 
.788699E 
.006000E 
.137054E 
.271038E 
.406972E 
.543897E 
.680882E 
.817027E 
.951472E 
.083396E 
.212024E 
.336632E 
.456543E 
.571135E 
.679838E 
.782139E 
.877577E 
.965748E 
.046304E 
.118949E 
.183441E 
.239591E 
.287260E 
.326358E 
.356843E 
.378717E 
.392028E 
.396861E 
.393345E 
.381640E 


06 
06 
06 
06 
06 
06 
06 
05 
05 
05 
05 
05 
05 
05 
05 
05 
04 
04 
04 
03 
02 
04 
04 
04 
05 
05 
05 
05 
05 
05 
05 
05 
05 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
06 
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1.162357E-08 

.833561E 

03 

3. 

.361946E 

06 

1.172357E-08 

.826059E 

03 

3. 

.334490E 

06 

1.182357E-08 

.816461E 

03 

3. 

.299532E 

06 

1.192357E-08 

.804815E 

03 

3. 

.257357E 

06 

1.202357E-08 

.791166E 

03 

3. 

.208274E 

06 

1.212357E-08 

.775561E 

03 

3< 

.152617E 

06 

1.222357E-08 

.758049E 

03 

3. 

.090738E 

06 

1.232357E-08 

.738680E 

03 

3. 

.023007E 

06 

1.242357E-08 

.717501E 

03 

2. 

.949809E 

06 

1.252357E-08 

.694563E 

03 

2. 

.871543E 

06 

1.262357E-08 

.669916E 

03 

2. 

.788618E 

06 

1.272357E-08 

.643610E 

03 

2. 

.701454E 

06 

1.282357E-08 

.615697E 

03 

2. 

.610476E 

06 

1.292357E-08 

.586227E 

03 

2. 

.516115E 

06 

1.302357E-08 

►555250E 

03 

2i 

.418804E 

06 

1.312357E-08 

.522819E 

03 

2. 

,318978E 

06 

1.322357E-08 

.488984E 

03 

2. 

.217074E 

06 

1.332357E-08 

.453796E 

03 

2. 

.113524E 

06 

1.342357E-08 

.417306E 

03 

2. 

.008757E 

06 

1.352356E-08 

.379565E 

03 

1. 

►903200E 

06 

1.362356E-08 

.340623E 

03 

1 

.797270E 

06 

1.372356E-08 

.300530E 

03 

1. 

.691379E 

06 

1.382356E-08 

.259337E 

03 

1. 

.585930E 

06 

1.392356E-08 

.217093E 

03 

1. 

►481316E 

06 

1.402356E-08 

.173848E 

03 

1. 

.377919E 

06 

1.^12356E-08 

.129651E 

03 

1. 

.276110E 

06 

1.^22356E-08 

.084549E 

03 

1 

.176247E 

06 

1.432356E-08 

.038592E 

03 

1. 

.078674E 

06 

1.442356E-08 

9. 

.918276E 

02 

9. 

.837220E 

05 

1.^52356E-08 

9. 

.443025E 

02 

8 

►917072E 

05 

1.462356E-08 

8 

.960637E 

02 

8 

.029302E 

05 

1.472356E-08 

8. 

.471576E 

02 

7. 

.176760E 

05 

1.482356E-08 

7 

.976300E 

02 

6 

.362136E 

05 

1.492356E-08 

7 

.475260E 

02 

5. 

.587952E 

05 

1.502356E-08 

6. 

►968905E 

02 

4 

.856564E 

05 

1.512356E-08 

6 

.457675E 

02 

4. 

.170156E 

05 

1.522356E-08 

5. 

.942005E 

02 

3. 

.530742E 

05 

1.532356E-08 

5. 

.422325E 

02 

2. 

.940161E 

05 

1.542356E-08 

4. 

.899058E 

02 

2. 

.400077E 

05 

1.552356E-08 

4< 

.372622E 

02 

1. 

.911982E 

05 

1.562356E-08 

3. 

.843427E 

02 

1. 

.477193E 

05 

1.572356E-08 

3. 

.311878E 

02 

1« 

096854E 

05 

1.582356E-08 

2. 

.778373E 

02 

7, 

.719356E 

04 

1.592356E-08 

2. 

.243303E 

02 

5< 

►032410E 

04 

1.602356E-08 

1. 

.707055E 

02 

2« 

.914035E 

04 

1.612356E-08 

1. 

.170005E 

02 

1. 

.368911E 

04 

1.622356E-08 

6. 

.325258E 

01 

4. 

.000888E 

03 

1.632356E-08 

9. 

.498233E 

00 

9< 

.021642E 

01 

1.642356E-08 

-4, 

.422673E 

01 

1. 

.956004E 

03 

1.652356E-08 

-9. 

.788719E 

01 

9, 

.581901E 

03 

1.662356E-08 

~1. 

.514487E 

02 

2. 

.293670E 

04 

1.672356E-08 

-2i 

.048775E 

02 

4. 

.197478E 

04 

1.682356E-08 

-2i 

.581405E 

02 

6< 

.663652E 

04 

1.692356E-08 

-3. 

.112054E 

02 

9. 

.664880E 

04 

1.702356E-08 

-3. 

.640405E 

02 

1< 
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9.842348E-03  2.589285E  03  6.704395E  06 

9.942348E-08  2.515993E  03  6.330221E  06 

1.004235E-07  2.443669E  03  5.971516E  06 

1.014235E-07  2.372375E  03  5.628165E  06 

1.024235E-07  2.302172E  03  5.299995E  06 

1.034235E-C7  2.233109E  03  4.986776E  06 

1.044235E-07  2.165232E  03  4.688230E-  06 

1.054235E-07  2.098581E  03  4.404043E  06 

1.064235E-07  2.033191E  03  4.133865E  06 

1.074235E-07  1.969090E  03  3.877316E  06 

1.084235E-07  1.906304E  03  3.633994E  06 

1.094235E-07  1.844851E  03  3.403476E  06 

1.104235E-07  1.784748E  03  3.185326E  06 

1.114235E-07  1.726006E  03  2.979098E  06 

1.124235E-07  1.668634E  03  2.784339E  06 

1.134235E-07  1.612635E  03  2.600593E  06 

1.144235E-07  1.558013E  03  2.427403E  06 

1.154235E-07  1.504764E  03  2.264315E  06 

1.164235E-07  1.452886E  03  2.110877E  06 

1.174235E-07  1.402371E  03  1.966645E  06 

1.184235E-07  1.353211E  03  1.831180E  06 

1.194235E-07  1.305395E  03  1.704057E  06 

1.204235E-07  1.258910E  03  1.584855E  06 

1.214235E-07  1.213742E  03  1.473170E  06 

1.224235E-07  1.169875E  03  1.368607E  06 

1.234235E-07  1#127290E  03  1.270784E  06 

1.244235E-07  1.085971E  03  1.179333E  06 

1.254235E-07  1.045897E  03  1.093900E  06 

1.264234E-07  1.007047E  03  1.014143E  06 

1.274234E-07  9.693999E  02  9,397362E  05 

1.284234E-07  9.329338E  02  8.703655E  05 

1.294234E-07  8.976253E  02  8.057312E  05 

1.304234E-07  8.634509E  02  7.455474E  05 

1.314234E-07  8.303867E  02  6.895421E  05 

1.324234E-07  7.984085E  02  6.374561E  05 

1.334234E-07  7.674917E  02  5.890435E  05 

1.344234E-07  7.376115E  02  5.440707E  05 

1.354234E-07  7.087430E  02  5.023166E  05 

1.364234E-07  6.808612E  02  4.635720E  05 

1.374234E-07  6.539412E  02  4.276391E  05 

1.384234E-07  6.279581E  02  3.943313E  05 

1.394234E-07  6.028870E  02  3.634727E  05 

1.404234E-07  5.787032E  02  3.348974E  05 

l,414234E-07  5.553823E  02  3.084495E  05 

1.424234E-07  5.329000E  02  2.839824E  05 

1.434234E-07  5.112324E  02  2.613686E  05 

1.444234E-07  4.903560E  02  2.404490E  05 

1.454234E-07  4.702472E  02  2.211325E  05 

1.464234E-07  4.508830E  02  2.032955E  05 

1.474234E-07  4.322407E  02  1.868320E  05 

1.484234E-07  4.142979E  02  1.716427E  05 

1.494234E-07  3.970327E  02  1.576349E  05 

1.504234E-07  3.804236E  02  1.447221E  05 

1.514234E-07  3.644495E  02  1.328234E  05 


1.52423AE-07 

1.534234E-07 

1.54A234E-07 

1.554234E-07 

1.564234E-07 

1.574234E-07 

1.584234E-07 

1.594234E-07 

1.604234E-07 

1.614234E-07 

1.624234E-07 

1.634234E-07 

1.644234E-07 

1.654234E-07 

1.664234E-07 

1.674234E-07 

1.684234E-07 

1.694234E-07 

1.704234E-07 

1.714234E-07 

1.724234E-07 

1.734234E-07 

1.744234E-07 

1.754234E-07 

1.764234E-07 

1.774234E-07 

1.784234E-07 

1.794234E-07 

1.804234E-07 

1.814234E-07 

1.824234E-07 

1.834234E-07 

1.844233E-07 

1.854233E-07 

1.864233E-07 

1.874233E-07 

1.884233E-07 

1.894233E-07 

1.904233E-07 

1.914233E-07 

1.924233E-07 

1.934233E-07 

1.944233E-07 

1.954233E-07 

1.964233E-07 

1.974233E-07 

1.984233E-07 

1.994233E-07 

2.004233E-07 

2.014233E-07 

2.024233E-07 

2.034233E-07 

2.044233E-07 

2.054233E-07 

2.064233E-07 

2.074233E-07 

2.084233E-07 

2.09^233E-07 

2.104233E-07 

2.114233E-07 

2.124233E-07 

2.134233E-07 
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3.490897E  02  - 
3.343241E  02 
3.201329E  02 
3.064969E  02 
2.933970E  02 
2.808150E  02 
2«687329E  02 
2.571332E  02 
2.459990E  02 
2.353138E  02 
2.250614E  02 
2.152263E  02 
2.057932E  02 
1.967473E  02 
1.880745E  02 
1.797607E  02 
1.717926E  02 
1.641571E  02 
1.558415E  02 
1.498337E  02 
1.431219E  02 
1.366946E  02 
1.305407E  02 
1.246493E  02 
1.190109E  02 
1.136147E  02 
1.084513E  02 
1.035115E  02 
9.878620E  01 
9.426686E  01 
8.994511E  01 
8.581291E  01 
8.186249E  01 
7.808642E  01 
7.447750E  01 
7.102883E  01 
6.773373E  01 
6.458580E  01 
6.157888E  01 
5.870707E  01 
5.596465E  01 
5.334617E  01 
5.084635E  01 
4.846013E  01 
4.618268E  01 
4.400932E  01 
4.193558E  01 
3.995717E  01 
3.806997E  01 
3.627004E  01 
3.455357E  01 
3.291694E  01 
3.135667E  01 
2.986942E  01 
2.845201E  01 
2.710137E  01 
2.581457E  01 
2.458880E  01 
2.342138E  01 
2.230973E  01 
2.125141E  01 
2.024404E  01 


1.218636E  05 
1.117726E  05 
1.024851E  05 
9.394034E  04 
8.608182E  04 
7.885707E  04 
7.221736E  04 
6.611747E  04 
6.051551E  04 
5.537257E  04 
5.065263E  04 
4.632234E  04 
4.235082E  04 
3.870952E  04 
3.537201E  04 
3.231391E  04 
2.951269E  04 
2.694754E  04 
2.459926E  04 
2.245015E  04 
2.048388E  04 
1.868540E  04 
1.704086E  04 
1.553751E  04 
1.416359E  04 
1.290829E  04 
1.176168E  04 
1.071462E  04 
9.758713E  03 
8.886242E  03 
8.090123E  03 
7.363855E  03 
6.701468E  03 
6.097489E  03 
5.546898E  03 
5.045095E  03 
4,587858E  03 
4.171325E  03 
3.791959E  03 
3.446520E  03 
3.132042E  03 
2.845814E  03 
2.585351E  03 
2.348384E  03 
2.132839E  03 
1.936820E  03 
1.758593E  03 
1.596575E  03 
1.449323E  03 
1.315516E  03 
1.193949E  03 
1.083525E  03 
9.832405E  02 
8.921825E  02 
8.095171E  02 
7.344845E  02 
6.663922E  02 
6.046093E  02 
5.485611E  02 
4.977242E  02 
4.516222E  02 
4.098212E  02 


.14423 
.15423 
.16423 
.17423 
.18423 
.19423 
.20423 
.21423 
.22423 
.23423 
.24423 
.25423 
.26423 
.27423 
.28^23 
.29423 
.30423 
.31423 
.32423 
.33423 
.34423 
.35423 
.36423 
.37423 
.38423 
.39423 
.40423 
.41423 
.42423 
.43423 
.44^23 


3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
3E-07 
2E-07 
2E-07 
2E-07 
2E-07 


1.928 
1.837 
1.750 
1.668 
1.589 
1.515 
1.444 
1.377 
1.313 
1.252 
1.195 
1.140 
1.088 
1.039 
9.935 
9.497 
9.083 
8.692 
e-323 
7.976 
7.b48 
7.341 
7.052 
6.781 
6.528 
6.292 
6.073 
5.869 
5.680 
5.507 
5.348 


539E 
331E 
574E 
073E 
640E 
097E 
273E 
006E 
141E 
530E 
033E 
515E 
849E 
915E 
966E 
849E 
765E 
730E 
815E 
140E 
878E 
245E 
504E 
961E 
963E 
904E 
214E 
364E 
861E 
253E 
121E 


26  . 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 


3.71926 
3.37578 
3.06450 
2.78246 
2.52695 
2.29551 
2.08592 
1.89614 
1.72434 
1.56883 
1.42810 
1.30077 
1.18559 
1.08142 
9.87234 
9.02091 
8.25147 
7.55635 
6.92858 
6.36188 
5.85053 
5.38938 
4.97378 
4.59949 
4.26273 
3.96006 
3.68839 
3,44494 
3.22721 
3.03298 
2.86024 


3E 
4E 
9E 
7E 
4E 
8E 
5E 
6E 
OE 
3E 
4E 
4E 
3E 
3E 
2E 
4E 
8E 
5E 
9E 
IE 
3E 
8E 
2E 
9E 
6E 
4E 
3E 
3E 
8E 
4E 
OE 


02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
01 
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